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1.  INTRODUCTION 

Let  Z^,  Z2,...,  Z^  be  a  random  sample  of  size  n  from  the  logistic  L(0,^)  population 
with  probability  density  function 

f  (z)  =  e~*/(l+e~*)^,  -ffl  <  z  <  »,  (1*1) 


-a  <  Z  <  », 


and  cumulative  distribution  function 
F*(z)  =  1/(1  +  e-*). 


-®  <  z  <  ®. 


Let  Zj,^  <  Z2.^  <...<  Z^.j^  denote  the  order  statistics  obtained  by  arranging  the  above 
sample  in  increasing  order  of  magnitude.  Then  the  density  function  of  Z..^  (l<i<n)  is 


given  by 


-®<Z.<B, 
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and  the  joint  density  function  of  Z-  and  Z.  (l<i<j<n)  is  given  by 

l.Il  J*ll  “  * 


¥  w 

{  (z.)  f  (z.),  -0D<Z.<Z.<B. 

«l  J 


(1.4) 


Let  us  now  denote  the  single  moments  E(Z^.^)  by  for  l<i<n,  k>l,  and  the  product 
moments  by  a*  for  l<i<j<n.  For  convenience,  let  us  also  denote 

by  a*,j^  and  E(Z?,^)  by  a*  for  l<i<n.  Further,  let  us  denote  Cov(Zj.jj^,Zj,j^)  by  0*^ 

Order  statistics  Z-.^  and  their  moments  have  been  studied  in  great  detail  by  several 
authors  including  Bimbaum  and  Dudman  (1963),  Gupta  and  Shah  (1965),  Tarter  and 
Clark  (1965),  Shah  (1966,  1970),  Gupta  et  al.  (1967),  Malik  (1980),  George  and  Rousseau 
(1987),  and  Balakrishnan  and  Malik  (1990).  Bimbaum  and  Dudman  (1963)  derived 
explicit  expression  for  the  cumulants  of  order  statistics  and  tabulated  the  means  and 
standard  deviations  for  sample  sizes  up  to  ten  and  for  some  large  sample  sizes  as  well. 
They  then  summarized  these  quantities  in  graphs  to  facilitate  interpolation  to  other 
sample  sizes.  Gupta  and  Shah  (1965)  derived  exact  expressions  for  the  moments  of  order 
statistics  in  terms  of  Bernoulli  and  Stirling  numbers  of  first  kind  and  used  them  to 
tabulate  the  first  four  moments  for  sample  sizes  up  to  ten.  They  also  expressed  the 
cumulants  in  terms  of  polygamma  functions,  as  was  originally  pointed  out  by  Plackett 
(1958).  It  should  be  mentioned  here  that  Plackett  (1958)  used  these  explicit  expressions 
of  the  moments  of  logistic  order  statistics  to  develop  a  method  of  approximating  the 
moments  of  order  statistics  from  an  arbitrary  continuous  distribution.  Distribution  of  the 
sample  range  has  been  studied  by  Gupta  and  Shah  (1965)  who  also  provided  a  short  table 
of  its  percentage  points  for  n=2  and  3.  By  generalizing  this  result,  Malik  (1980)  derived 
the  exact  formula  for  the  cumulative  distribution  function  of  the  r^^  quasi-^ange,  viz.. 
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^n-r  n'’^r+l  n  ’  [^] '  “  independent  study,  Tarter  and  Clark  (1965) 

reproduced  some  of  the  results  of  Gupta  and  Shah  (1965)  and  then  studied  the 
distribution  of  the  sample  median  in  detail.  George  and  Rousseau  (1987)  recently 
examined  the  distribution  of  the  sample  midrange,  viz.,  (Zi._+Z  )/2,  and  established 

several  relationships  in  distribution  between  the  midrange  and  sample  median  of  the 
logistic  and  Laplace  random  variables.  A  series  expression  for  the  covariance  of  two  order 
statistics  has  been  provided  by  both  Gupta  and  Shah  (1965)  and  Tarter  and  Clark  (1965). 
Shah  (1966)  tabulated  the  covariances  for  sample  sizes  up  to  ten  and  Gupta,  Qureishi  and 
Shah  (1967)  extended  this  table  for  sample  sizes  up  to  twenty  five.  It  should  be 
mentioned  here,  however,  that  by  means  of  some  recurrence  formulas  Kjelsberg  (1962) 
had  already  derived  exact  numerical  results  for  the  covariances  jErom  samples  of  size  five 
or  less. 

By  using  the  fact  that  f  (z)  and  F*(z)  given  in  Eqs.  (1.1)  and  (1.2)  satisfy  the 
rdation 


t’(.)  =  P*w{l-F*(*)}.  (1.5) 

Shah  (1966,  1970)  established  several  recurrence  relations  satisfied  by  the  single  and  the 
product  moments  of  order  statistics.  Recently,  Balakrishnan  and  Vilik  (1990)  prepared 
tables  of  means,  variances  and  covariances  for  sample  sizes  up  to  i  r  oy  applying  these 
rdations  in  a  simple  and  systematic  recursive  way.  Some  of  the  results  in  the  references 
dted  above  have  also  been  summarized  in  a  review  article  by  Malik  (1985). 

The  truncated  logistic  distribution  plays  a  role  in  a  variety  of  applications,  as  has 
been  mentioned  by  Kjelsberg  (1962).  Order  statistics  and  their  moments  from  a  general 
truncate  logistic  distribution  have  been  studied  by  Tarter  (1966).  He  derived  exact  and 
explicit  expressions  for  the  means,  variances  and  covariances  of  order  statistics  in  terms 
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of  a  finite  series  involving  logarithms  and  dilogarithms  of  the  constants  of  truncation.  By 
following  the  lines  of  Shah  (1966,  1970),  Balakrishnan  and  Joshi  (1983  a,b)  established 
several  recurrence  relations  satisfied  by  the  single  and  the  product  moments  of  order 
statistics  from  a  symmetrically  truncated  logistic  distribution.  Then,  Balakrishnan  and 
Kocherlakota  (1986)  generalized  these  results  to  the  doubly  truncated  logistic  distribution 
and  displayed  that  these  recurrence  relations  could  be  used  systematically  in  order  to 
evaluate  the  means,  variances  and  covariance  of  all  order  statistics  for  all  sample  sizes. 

In  this  paper  we  present  a  detailed  discussion  of  order  statistics  from  the  logistic 
distribution  and  some  of  their  properties.  In  Section  2  we  give  the  percentage  points  and 
modes  of  order  statistics.  In  Section  3  we  derive  exact  and  explicit  expressions  for  the 
single  and  the  product  moments  of  order  statistics.  These  work  in  terms  of  gamma 
function  and  its  successive  derivatives.  In  Section  4  we  present  some  recurrence  relations 
satisfied  by  the  single  and  the  product  moments  of  order  statistics  which  would  enable 
one  to  compute  the  means,  variances  and  covariances  in  a  simple  recursive  way.  The 
distribution  function  of  the  sample  range,  as  derived  by  Gupta  and  Shah  (1965),  is 
presented  in  Section  5  and  the  distribution  of  the  quasi^ange  derived  by  Malik 
(1980)  is  also  given  here  for  the  sake  of  completeness.  In  Section  6  we  give  some  relations 
between  moments  for  the  case  of  the  doubly  truncated  logistic  distribution  that  are  due 
to  Balakrishnan  and  Kocherlakota  (1986).  In  Section  7  we  present  details  of  tables  that 
are  available  in  this  context.  Finally,  in  Section  8  we  describe  Blackett’s  (1958)  method 
of  approximating  the  moments  of  order  statistics  from  an  arbitrary  continuous 
distribution  by  using  the  moments  of  logistic  order  statistics. 

2.  PERCENTAGE  POINTS  AND  MODES 

The  distribution  function  of  is  given  by 


(2.1) 


From  (2.1)  we  obtain  the  100a  percentage  point  of  to  be 


.  0<a<l.  (2.2) 

Next,  the  distribution  function  of  Z  is  given  by 

Xa*U 

CW)  =  {''*(*«)}“  (2-3) 

From  (2.3)  we  obtain  the  100a  percentage  point  of  Z^.^  to  be 

=  0<o<l.  (2.4) 

Similarly,  the  distribution  function  of  (2<i<n-l)  is  given  by 

=  B(i,ii4+1)  J  “  (M  (2-5) 

0 

Now,  let  B^(i,n-4+l)  denote  the  100a  percentage  point  of  the  Beta(i,n-4+l)  distribution. 
From  (2.5)  we  then  have  the  100a  percentage  point  of  Zj.^^  (2<i<n-l)  to  be 

*i:n(a)  =  ^  Ji,n-i4-l)},  0<a<l.  (2.6) 
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By  using  the  symmetric  relation  satisfied  by  the  incomplete  beta  functions,  we  observe 
firom  Eqs.  (2.2),  (2.4)  and  (2.6)  that  the  100a  percentage  point  of  Z..^  is  simply  the 
negative  of  lOO(l-o)  percentage  point  of  While  the  percentage  points  of 

and  Z^.^  may  be  obtained  easily  from  (2.2)  and  (2.4),  respectively,  the  percentage  points 
of  Z|.^  for  2<i<n-l  may  be  obtained  from  (2.6)  either  by  using  the  extensive  tables  of 
incomplete  beta  function  prepared  by  Karl  Pearson  (1934)  and  Pearson  and  Hartley 
(1970)  or  by  using  the  algorithm  given  by  Cran,  Martin  and  Thomas  (1977).  Gupta  and 
Shah  (1965)  have  tabulated  some  percentage  points  of  all  order  statistics  for  sample  sizes 
up  to  10  and  some  selected  order  statistics  for  sample  sizes  up  to  25. 

Next,  by  differentiating  the  density  function  of  Z..^  in  (1.3)  with  respect  to  Zj  and 
using  the  relation  in  (1.5),  we  get 


df.  (z.) 
1  :n'  r 


■HiT 


=  (i:7)T(n^T  ^*i)}  (*i)}  [Kn+l)^  (^i)]  • 


(2.7) 


Upon  equating  (2.7)  to  zero  and  solving  for  z^,  we  obtain  the  mode  of  Z-.^  (l<i<n)  to  be 

m*:n  =  Ai|i/(n-i+l)J  =  (a  (pj/q^)  ,  (2.8) 

where  Pj=l-qj=i/(n+l).  Due  to  the  symmetry  of  the  logistic  distribution,  we  observe 
once  again  that  the  mode  of  Z..^  is  simply  negative  of  the  mode  of  ^n-i+l:n' 
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3.  MOMENTS  AND  CUMULANTS 


From  (1.3)  we  have  the  moment  generating  function  of  (l<i<n)  to  be 

1  f  g-(n-i+l)z+t2 

=  b(i,n-i+l)  J  (i!j.e“®)n+l  ^ 


=  B(i+t,n-i+l^)/B(i,n-i+l) 


r(i+t)  r(n— i+1— t) 
-ifTT  ¥1)  ’ 


where  B(.,.)  and  r(.)  are  the  usual  complete  beta  and  gamma  functions,  respectively.  An 
alternate  expression  of  M- .^(t)  involving  Bernoulli  numbers  and  Stirling  numbers  of  first 
kind  has  been  given  by  Gupta  and  Shah  (1965).  From  the  expression  of  the  moment 
generating  function  in  (3.1),  we  obtain  the  following: 

"  ^^+1),  (3.2) 

“im  =  +  (iKi)  -  .  (3.3) 


where 


^i,i:n  =  +  ^'(““*+1)  . 

^z)  =  ^  Air(z)  = 


and 
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are  the  digamma  and  trigamma  functions,  respectively.  Thus,  from  Eqs.  (3.2)  -  (3.4), 
one  may  compute  the  means  and  variances  of  order  statistics  either  by  using  the 
extensive  tables  of  digamma  and  trigamma  functions  prepared  by  Davis  (1935)  and 
Abramowitz  and  Stegun  (1965)  or  by  using  the  algorithms  given  by  Bernardo  (1976)  and 
Schneider  (1978).  Gupta  and  Shah  (1965)  have  given  exact  expressions  for  the  first  four 
moments  of  order  statistics  for  sample  sizes  up  to  10  and  the  values  of  mean  and  variance 
have  been  tabulated  recently  by  Balakrishnan  and  Malik  (1990)  for  sample  sizes  up  to  50. 
We  may  note  here  that  the  moment  generating  function  in  (3.1)  may  be  used  to  obtain 
higher  order  single  moments  also  by  involving  polygamma  functions. 

From  (3.1)  we  get  the  cumulant  generating  function  of  (l<i<u)  to  be 


a=  Ai  r(i+t)  -f  Ai  r(n-i+l-4)  -  Ai  r(i)  -Ai  r(n-i+l).  (3.5) 


.th 


From  (3.5)  we  obtain  the  k  cumulant  of  (l<i<n)  to  be 


k  k 

=  At  ^  r(i+t)  I  +  Ai  r(n-i+l-t)  I 

dt^  U=0  dt* 


=  +  (-1)^  t^(^"^)(n-i+l), 


(3.6) 


where 


and 


.(k-1) 


(z)  =^Air(z) 
dz 


i^%)  =  Hz)  =  . 


for  k=l,2,.... 


It  is  clear  from  (3.6)  that  for  k=l,2,..., 

*(2k-l)  _  *(2k-l) 

^i:n  \-i+l:n 

(3,7) 

and 

^♦(2k)  _  ^♦(2k) 

''i:n  '*n-i+l:n‘ 

(3.8) 

These  may  also  be  observed  simply  by  using  the  symmetry  of  the  logistic  distribution. 
By  applying  the  series  expansions  for  ^(z)  and  given  by 


^=1 

and 

i^-%)  =  (-D^k-l).  I  -1  .  k>2. 

1=1 

we  obtain  from  (3.6)  that  for  n-i+l>i 

(3.9) 

(3.10) 

The  above  formulae  for  the  first  four  cumulants  were  originally  given  by  Plackett  (1958). 
From  (3.9)  we  get 


*l:a  =-{l+!  +  J+-+5^} 
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which  was  also  given  by  Gumbd  (1958).  Tbe  exact  and  explicit  expression  for  the 
cumulants  of  logistic  order  statistics  given  in  Eqs.  (3.9)  and  (3.10)  will  later  be  used  in 
Section  8  for  developing  some  series  approximations  for  the  moments  of  order  statistics 
from  an  arbitrary  continuous  distribution. 

From  the  joint  density  of  and  (l<i<j<n)  in  (1.4),  we  have  the  joint 
moment  generating  function  of  Z-.^  and  Zj.^  to  be 


M 


±  f  tiZ*  \ 


X  Z 


H-1  r  ♦.  .iH 


*  * 


f  (zj)  f  (Zj)  dZj  dzj. 


(3.11) 


Making  the  transformations 


u=F*(zj)  = - ^  and 

l+e  * 


v=F*(z.) 


and  thence,  noting  that 


and 
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we  can  rewrite  Eq.  (3.11)  as 


n! 


V  t,  t, 

f  u  ^ 

.  J  t,  t, 

0  0  (1-u)  ^  (1-v) 


u*  ^  ^  du  dv. 


-t 


’'1 

By  expanding  (1-u)  as  an  infinite  series  in  powers  of  u,  we  obtain 


M*  ft  t  1  y  (ti+M)^^  r 

^i.j:n(V2)  -  (i-l)!(j-i-l)!(n-i}!  Z  - Z -  J 


^=0 


0  0 


t,+i-H-/  j:  i  t« 
u  ^  (v-u)^  ^  V  ^(1-v)  ^  du  dv, 


where 


(t^+/-l)(^ 


=  1 


if  ^=0 


=  t  j(tj+l) (tj+M)  if  ^1. 


By  noting  that 


T  t,+i-l+/  1  j+t,+/-l 

J  u  ^  (v-u)^  ^  du  =  V  B(tj+i+/,}-i), 

0 


(3.12) 


we  may  rewrite  Eq.  (3.12)  as 
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=  (i-T)!(>-l-l)!(ii-j)!  I  - Z! - 

i=o 


0 


j+ti+tn+Z-l  n— j-t,j 

V  ^  2  2^^ 


;  (ti+«)W 

(i-l)!(H-l)!(ii-j)!  2.  - 2 - B(tj+iHj-i) 

Z=0 


B(  j+t  j^+t2+Z,n— j-t2+ 1) 

r(n+i )  ?  r(tj+i+z) 

r(i)r(n--j4-i)  Z  — 2 — r(t,+j+z) 
z=o  ^ 


r(tj+t2+j+/)r(n-j+i-t2) 

l'(n+l+tj+Z) 


(3.13) 


From  the  above  expression  for  the  joint  moment  generating  function  of  Z..^  and  Zj 
can  obtain  the  product  moments  as  follows: 


ki  k2 

o.  ^  =  E(Z.  Z.  )  = 
i,j:n  '  i:n  j:n' 


at, 


dt 


lE^^i,j:n(h’*2) 


ti-t2-0 


,  one 


(3.14) 


The  case  k^=k2=l  is  of  particular  importance  and,  in  this  case,  we  get 
“i,j:n  =  ^  ■  V(n+1)}  |v(j)  -  V<H+1)} 

^  ~  V(n-j+l)}  .  (3.15) 

Shah  (1966)  tabulated  the  covariances  of  order  statistics  for  sample  sizes  up  to  ten 
while  Gupta,  Qureishi  and  Shah  (1967)  extended  up  to  25.  Recently,  Balakrishnan  and 
Malik  (1990)  provided  tables  of  means,  variances  and  covariances  for  sample  sizes  50  and 
less.  It  should  also  be  mentioned  here  that  Balakrishnan  and  Leung  (1988  a,b)  derived 
series  expressions  similar  to  the  ones  given  in  Eqs.  (3.2)-(3.4)  and  (3.15)  for  the  single 
and  the  product  moments  of  order  statistics  from  a  generalized  logistic  distribution  and 
provided  tables  of  means,  variances  and  covariances  for  sample  sizes  up  to  15. 

4.  RECURRENCE  RELATIONS  FOR  MOMENTS 

In  this  section  we  shall  present  some  recurrence  relations  satisfied  by  the  single  and 
the  product  moments  that  were  established  by  Shah  (1966,  1970)  and  show  that  one  may 
evaluate  these  moments  for  all  order  statistics  from  all  sample  sizes  in  a  simple  and 
systematic  recursive  marjici. 


(4.1) 
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Proof.  From  (1.3)  we  have  for  n>l  and  k>0 


“im 


(D 


=  n  I  |i-F*(2)|^  ^  f* (z)  dz 


By  using  the  relation  (1.5)  in  Eq.  (4.2),  we  get 


a 


l:n 


QD  ^ 

=  n  I  z^  F  (z)  |l-F  (z)|  dz 


which,  upon  integrating  by  parts,  yields 


OD 


=  E7T  [“  j  ***(*)  {l-F*(z)|  f*(z)  dz 


-  I  |l-F*(z)j”  f* (z)  dz 


The  recurrence  relation  in  (4.1)  is  obtained  by  rewriting  the  above  equation. 


Relation  4.2:  For  l<i<n  and  k=0,l,2, . , 


“i+l:n+l  “l;n+l  ^  i  (  n-l+1 )  “i:n 


(4.2) 


(4.3) 
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Proof.  From  (1.3)  we  have  for  l<i<n  and  k>0 


“i^n'  =  J  {*'  W}  * 


(4.4) 


By  using  the  relation  (1.5)  in  Eq.  (4.4),  we  get 


which,  upon  integrating  by  parts,  yields 


=  (i-l)!(n-i)!(k+l)  [(“-*+!)  j  {f  (*)}  (l-*'  (2)}  ^  (*)d* 


-i}zl^^l{F%)f-^(l-F%)pS?z)dz] 


_  i(n-i+l)  f.*(k+l)  .♦(k+l)l 

-  (k+l)(n+l)  t®i+l :  n+1  "  ®i:n+l  / 


The  recurrence  relation  in  (4.3)  follows  by  rewriting  the  above  equation. 

With  the  values  of  (j=l,2,...,k)  known,  one  may  be  able  to  use  Relations  4.1 

and  4.2  in  a  simple  recursive  way  to  compute  the  first  k  single  moments  of  all  order 

statistics  from  all  sample  sizes.  Thus,  for  example,  by  starting  with  the  values  of 
’•'fll  *(2^  2 

and  employ  Relations  4.1  and  4.2  to  evaluate  the  first 
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two  single  moments  and,  thence,  the  variances  of  all  order  statistics  from  all  sample  sizes 
in  a  simple  and  systematic  recursive  process.  These  computations  may  be  checked  by 
using  the  identities  (David,  1981,  p.  39;  Arnold  and  Balakrishnan,  1989,  p.  6) 


n 


i=l 

see  also  Balakrishnan  and  Malik  (1986). 


(4.6) 


Relation  4.3:  For  l<i<n-l, 


a*  =  aPl  +  -SiiJo*  o*  1  (4.6) 

i,i+l:n+l  i:n+l  n-i  +1)  i,i+l:n  i:n  n-^  i:nj  .  ' 


Proof.  For  l<i<n-l,  we  may  write  from  (1.4)  that 


“i:n  =  J 


=  (i-l)!(n-i-l)!  1  *1  {'■  ^  (h)  ‘^1’ 


where 


H^l  1  n-i-1  4, 

l-F  (Zj)}  f  (zj)  dz 


(4.8) 
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By  using  the  relation  (1.5)  in  Eq.  (4.8)  and  integrating  by  parts,  we  get 

f  *  f  ♦  in-i+l 

K(zj)  =  -  Zj|l-F  (z^)|  +  Zj|l-F  (zj)j 

(P  1 

+  (n-i)  j  Zj  |i-F*(z2)}  f*(z2)  dz2 

*1 

(B 

-  (n-i+l)  I  Z2  {i-F*(z2)}  f* (Z2)  dz^  . 

*1 

Upon  substituting  the  above  expression  of  K(z^)  in  (4.7)  and  simplifying  the  resulting 
equation,  we  get 


The  recurrence  relation  in  (4.6)  follows  by  rewriting  the  above  equation. 
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Proof.  For  we  may  write  from  (1.4)  that 


a.  . ,  =  E(ZV  Z. . ,  ) 

i+l:n  '  i:n  i+l:n^ 


f* (2j)  K(*j)  <bj, 


li-1  Hn-i-D!  I  *2 


where 


*2 

K(*2)  =  J  d2j 


(4,10) 


(4.11) 


By  using  the  relation  (1.5)  in  Eq.  (4.11)  and  integrating  by  parts,  we  get 


K(*2)  =  ■*2{^*(*2)} 


^2  i— 1 

-i  j  *1  {f*(*i)} 


®2  I 

+  0+1)  j  *1  ^*(*i)  ^1 


Upon  substituting  the  above  expression  of  K(z2)  in  (4.10)  and  simplifying  the  resulting 
equation,  we  get 


i+l:n 


i«*(2)  i(i+l)  *(2) 

'  ®i+l:n  "  (n+1)  ®i+2:n+l 


,  ,  i(  i  4-1)  _* 

®i4+l:n  (n+l )  ®i+l,i+2;n+l’ 
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The  recuirence  relation  in  (4.9)  is  obtained  by  rewriting  the  above  equation. 
In  particular,  by  setting  i=n-l  in  Rdation  4.4  we  get  for  n>2 


a*  =a*(2)  .n+1 

“n,n+l:n+l  “n4-l:n+l  ^  n 


«  +  Vl,n:n  - 


It  should  be  mentioned  here  that  the  recurrence  relations  in  (4.6)  and  (4.12)  are 

sufficient  for  the  evaluation  of  all  the  product  moments  of  order  statistics  from  all  sample 

♦  *2 

sizes.  By  starting  with  the  result  that  2‘2~^1‘1~^  (Govindarajulu,  1963;  Joshi,  1971), 

the  recurrence  relations  in  (4.6)  and  (4.12)  will  enable  one  to  compute  all  the  immediate 

* 

upper-diagonal  product  moments  (l<i<n-I)  for  all  sample  sizes  in  a  simple 

* 

recursive  way.  All  the  remaining  product  moments,  viz.,  for  l<i<j<n  and  j-i>2, 
may  be  determined  systematically  by  employing  the  well-known  recurrence  relation 
(David,  1981,  p.  48;  Arnold  and  Balakrishnan,  1989,  p.lO) 


(4.13) 


that  is  true  for  any  arbitrary  distribution.  These  computations  may  then  be  checked  by 
using  the  identity  (David,  1981,  p.  39;  Arnold  and  Balakrishnan,  1989,  p.  10) 


I  I  ®i,j:n  ~ 

i=l  j=i+l 


(4.14) 


see  also  Balakrishnan  aT>d  Malik  (1986). 


By  proceeding  on  similar  lines,  we  may  also  establish  the  following  recurrence 
relations. 


Relation  4.5:  For  l<i<j<n  and  i-i>2, 


_♦  _ ♦  .  n+1  f.*  1  1 

®i,j:n+l  ”  ®i  j-l;n+l  n-j+2(®i,j:n  “  ®i,j-l;n  ”  n-j+1  ®i:nj' 


(4.15) 


Relation  4.6:  For  l<i<j<n  and  i-i>2, 


* 

“i+l,j+l:n+l 


* 

®i+2,j+l:n+l 


+ 


n+lfl 

mir 


+ 


M+i 


.j:4’ 


(4.16) 


It  should  be  pointed  out  here  that  one  may  employ  Relations  4.5  and  4.6  to 

determine  all  the  product  moments  other  than  the  immediate  upper-diagonal  product 
♦ 

moments,  viz.,  a.  for  l<i<j<n  and  instead  of  the  recurrence  relation  in  (4.13). 


5.  DISTRIBUTIONS  OF  SOME  SYSTEMATIC  STATISTICS 

In  this  section  we  first  present  the  distribution  of  the  sample  range  as  derived  by 
Gupta  and  Shah  (1965).  Then,  we  give  the  expression  of  the  distribution  of  the  r^^ 
quasi-range  derived  by  Malik  (1980). 

Let  us  denote  the  sample  range  Z  -  Z.  by  W  .  The  cumulative  distribution 

U*a&  X»tt  U 

function  of  can  be  written  down  as  (David,  1981,  p.  12) 


Pr(Wj^<w)  =  n  I  |f  (z+w)  -  F  (z)|  f  (z)  dz, 


0<w<(». 


(5.1) 
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{*  *  A  n-1 

F  (z+w)  -  F  (z)|  binomially,  we  get  from  (5.1)  that 


n— 1  _ 

Pt(W^<w)  =  n  J  (-1)’'  (”-*)  I  {f%+w)}’‘  {f’w}'  {%)  dz 

k=0 


OD 


n-l-k  f  ♦  ,  -jk  * 


n— 1 


-  “  I  (-1)“  (V)  I  (i+e-e-'r*-*  (i+e-^)»+^ 


k=0 


By  substituting  u=l/(l+e  ^e“*)  in  the  integral  in  (5.2),  we  get 


n— 1 


P'(W„<w)  =  B  J  (-!)'■  (Y)  6-0'+*)’'  Jw).  0<w<.,  (5.3) 

k=0 


where,  with  a  =  e  ^-1, 


du 


^  [(-l)'‘+l  (J7l)  *>(>+»)  +1^  (-!)'("?) 

/^k+1 


M 


[h)'  (i;})  -  A  ‘  (-1)'  (-f)  ■ 

0 

i^k+1 


(l-e-*) 


(5.4) 
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In  the  expression  of  Aj^  ^^(w)  in  (5.4),  should  be  set  to  zero  if  ^  *n-2.  By 
substituting  the  expression  of  A.  (w)  in  (5.4)  into  Eq.  (5.3),  we  derive  the  cumulative 

jE,n 

distribution  function  of  the  sample  range  as 

(V)[H)'‘  (kTl)  - 

^  k=0 

(-1)'  (V) 

/=0 

f#k+l 

(5.5) 


In  particular,  for  n=:2  and  3,  we  obtain  from  (5.5)  that 

Pr(W2<w)  =  |l  -  e”^^  -  2we'“^j/(l-e‘^)^  (5.6) 

and 

Pr(W3<w)  =  {l+9e"^-9e~^^-e"^^-6we"^(l+e"^)|/(l-e"^)^.  (5.7) 

Gupta  and  Shah  (1965)  tabulated  the  probability  integrals  of  the  range  for  n=2  and  3 
from  (5.6)  and  (5.7),  repsectively. 

By  differentiating  the  distribution  function  of  the  sample  range  in  (5.5)  with  respect 
to  w,  we  derive  the  density  function  of  the  sample  range  as 
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^y/  (^) - 


I  (-1)“  (V)  (kTl) 

f  k=0 


n— 1 


/=0 

/#k+l 


n— 1 


I  H)'‘  (”j‘)[(-l)'‘  (k7l)  e-<'‘+‘)''{l-{k+l)w} 


1=0 

/^k+1 

(5.8) 

where,  as  before,  should  be  set  to  zero  if  k>n-2. 

Proceeding  exactly  on  similar  lines,  Malik  (1980)  derived  the  cumulative 
distribution  function  of  the  r*^  quasi-range  r~^n-r-n  ”  ^r+l  n 


Pr(W^^^  <  w) 


2r-k 
r  n  (n-i) 
i=0 

FTTFITr 
k=0 


=  I 


n-2r +k-l 

J  (_l)j  |n-2r+k-lj  g-(r+j+l)w 
j=0  ^ 


r — k  /  |\/+1  fr — k 
I ^  ^J_g-wjn-r+k-l-^  '  [  r+j+1 


n— r  +  k+  f-1 
m=0 

m#  r  +  j+1 


g-(m-r-j-l)w_j 
m— r-j-1 


I 


0<W<OD, 


(5.9) 


-24- 


where  r+j+i  should  be  set  to  zero  if  j>n-2r+k+/-2.  The  distribution  function  of 
the  sample  range  in  (5.8)  may  be  derived  as  a  special  case  from  (5.9)  by  setting  r=0. 


6.  RESULTS  FOR  TRUNCATED  DISTRIBUTIONS 

In  this  section  we  start  with  order  statistics  from  a  doubly  truncated  logistic 

distribution  and  present  the  results  of  Balakrishnan  and  Joshi  (1983  a)  and  Balakrishnan 

and  Kocherlakota  (1986).  In  addition  to  generalizing  the  relations  given  in  Section  4, 

these  results  will  enable  one  to  evaluate  the  single  and  the  product  moments  of  all  order 

statistics  from  all  sample  sizes  in  a  simple  recursive  manner. 

Let  Z.  „<Z«  „< . <Z„.„  be  order  statistics  from  a  random  sample  of  size  n  from  a 

l.n~  2:n~  n:n  ^ 

doubly  truncated  logistic  distribution  with  probability  density  function 


where  Q  and  1-P  are  the  proportions  of  truncation  on  the  left  and  the  right  of  the 
standard  logistic  density  function  in  (1.1).  Under  this  notation, 


Qi=^[t^]« 


(6.3) 
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and  let 


From  (6.1)  and  (6.2)  we  observe  the  relations 

f**(2)  =  (l-2Q)F**(z)  -  (P-Q){P**W}=*  +  Qj.  (6.5) 

f*%)  =  (2P-1){1-P*%)}  -  (P-Q){1-p’*(2)}2  +  Pj,  (6.6) 

and 

f*%)  =  (2P-l)P**(z){l-F**(z))  +  (P+C!-1)(1-F**(z)}2  +  Pj. 

(6.7) 

Let  us  now  denote  the  single  moments  E(Zj  .^j)  by  Oj  '  (l<i<n,k>l)  and  the  product 

moments  E(Zj,^Zj.^)  by  Oj  (l<i<j<n).  For  convenience,  let  us  also  use  for 
**fll  **  **(2\ 

®i  ‘n  ®i  i  n  ®i  ‘n  '  Then,  these  moments  satisfy  the  following  recurrence 
relations. 


Relation  6.1:  For  k=0,l,2,.... 


+  pij[p4pr‘  -  q1[«)  +  (2P-l){a;*('‘«)- 


(6.8) 


J 
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Proof.  For  k>0,  let  us  consider 


^1" 


1  1.  ** 


=  z  f  (z)  dz 


Upon  using  (6.6)  in  the  above  equation  and  then  integrating  by  parts,  we  get 


pj(pJ+i-q;+'}]  . 


(6.9) 


The  recunence  relation  in  (6.8)  follows  by  rewriting  Eq.  (6.9). 


Relation  6.2:  For  k=0,l,2,...., 


-(i+i).;*!-)] . 


(6.10) 


Proof.  For  k>0,  let  us  consider 


1  k  ♦♦ 

=  z*  f  (z)  dz 
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Upon  using  (6.5)  in  the  above  equation  and  then  integrating  by  parts,  we  get 


“l:l 


-  (p-<5){p{+i-  •  (6.11) 


The  recurrence  relation  in  (6.10)  is  obtained  by  rewriting  Eq.  (6.11). 


J:  For  n>2  and  k=0,l,2,..., 


(6.12) 


Proof.  For  n>2  and  k>0,  let  us  consider 


Pi  1 

z*  |l  -  F  (z)|  f  (z)  dz 


Upon  using  (6.6)  in  the  above  equation  and  then  integrating  by  parts,  we  get 
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The  recurrence  relation  in  (6.12)  follows  by  rewriting  Eq.  (6.13). 


Relation  6.4:  For  k=0,l,2 . 


“7“  “2 


T] 


(6.14) 


Pmof.  For  k>0,  let  U8  consider 


“2:2 


'(■ 


2  z*  F  (z)  f  (z)  dz 


Upon  using  (6.6)  in  the  above  equation  and  then  integrating  by  parts,  we  get 


2Pj{p^+l  - 


The  recurrence  relation  in  (6.14)  is  obtained  by  rewriting  Eq.  (6.15). 


Relation  6.5:  For  n>3  and  k=0,l,2,...., 


•*(1+1)  =  o**(k+l)  ,  5±1  r!2J  "(l+l)  _  o**(k+l))  .  2Pzl 
“2:n+l  “l:n+l  +  F^[n^“2:n-1  “l:n-l  n 


(»(k+i).„;-(k+i)}.^^.(k)]. 


Proof.  For  n>3  and  k>0,  let  us  consider 
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Proof.  For  2<i<n-l  and  k>0,  let  us  consider 


.1®“*  ♦♦ 

f  (z)  dz 


Upon  using  (6.7)  in  the  above  equation  and  then  integrating  by  parts,  we  get 


a 


**W  ^  1 

i  :n  k+1 


(6.19) 


If  we  now  use  the  well-known  relation  (David,  1981,  p.  46;  Arnold  and  Balakrishnan, 
1989,  p.  6) 

- «!} = 

in  (6.19)  and  simplify  the  resulting  equation,  we  derive  the  recurrence  relation  in  (6.18). 


Relation  6.7:  For  n>2  and  k=0,l,2,..... 


♦♦(k+1)  _  n+1 


'n+l:n+l 


-5^(i>+1)(P+Q-1)-(2P-1)} 

+  (P+Q-1) • 


(6.20) 
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Proof.  For  n>2  and  k>0,  let  us  consider 


=  ■‘I''’'  {»'"(«)}"■*  '**«  ^ 


Upon  using  (6.7)  in  the  above  equation  and  then  integrating  by  parts,  we  get 


♦♦(k)  =  1  r  n  .2P_1^  L**(k+1)  _  «**(k+l)l 

“n:n  k+1  [n+1  '  ^pn+l:n+l  “n:n+l  J 


(6.21) 


If  we  now  use  the  well-known  relation  (David,  1981,  p.  46;  Arnold  and  Balakrishnan, 
1989,  p.  6) 


in  (6.21)  and  simplify  the  resulting  equation,  we  derive  the  recurrence  relation  in  (6.20). 

By  starting  with  the  values  of  (j=l>2,....,k),  one  can  employ  Relations 
6.1  -  6.7  in  a  simple  and  systematic  recursive  way  to  compute  the  first  k  single  moments 
of  all  order  statistics  from  all  sample  sizes  for  any  choice  of  Q  and  P. 

We  may  also  establish  the  following  recurrence  relations  satisfied  by  the  product 


moments  of  order  statistics. 
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Itelation  6.8:  We  have 


Proof.  Let  us  start  with 


♦♦(1)  _  Erz*  z*o  ^ 

"l:2  -^^h:2h:2f 


V  • 


(zj^)  K(zj^)  dzp 


where 


Pj 

K(*i)  “  J  ^  ^*2^  *^*2  ' 


(6.23) 


(6.24) 


Upon  using  (6.6)  in  (6.24)  and  integrating  by  parts,  we  get 


K(zj)  =  (2P-1)  j  Zj  Mzj)  ^2  ^  *1  {^  "  ^**(*l)} 


-  (P-<3)  2  j  Zj  |i-F**(z2)}  ^**(*2)  ^2  ~ 
•  «- 


+  P2  (^1  -  ‘1) 
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Upon  substituting  the  above  expression  of  K(*j)  in  (6.23)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.22). 


Relation  6.9:  We  have 


-(l-2<}){<4*|2)-a”4].  (6.25) 


Proof.  Let  us  start  with 


“2:2  ~®(^1:2^2:2^ 


'C 


=  2  [  *2  ^**(*2^  ^(*2^  *^2’ 


where 


K(*2)  =*  I  f**(*i)  d*l  ' 

Upon  using  (6.5)  in  (6.27)  and  integrating  by  parts,  we  get 

z  2 

K(*2)  =  (»-2Q){»2  *'**(*2)  -  ^  ‘1  '**1} 


(6.26) 


(6.27) 


(P-Q)[«2{F’*(«2)f  -  *'**(*1)  "(*l)  '‘*1] 


+  Q2 
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Upon  substituting  the  above  expression  of  K(z2)  in  (6.26)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.25). 


Relation  6.10:  For  l<i<n-2, 


O: 


i  ,i+l:n+l 


Proof.  For  l<i<n-2,  let  us  consider 


** 


(1) 


1  :n 


where 


=  *1  {*■  (‘l)}  *  ^1’ 


fl  r  **  ^n-i-1 

K(aj)  =  J  (*2)}  ^  *^2  ' 


(6.29) 


(6.30) 
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Upon  using  (6.6)  in  (6.30)  and  integrating  by  parts,  we  get 


K(*j)  =  (2P-l)[(ii-i)j\  f“(zj)  dZj  -Zj  {l-P**(Zj)}‘‘'^] 


-  (P-Q)[(n-i+l)]^j  {l-P*%2)}“‘’  dzj  -  {l-p‘*(zj)}'‘‘*'^‘ 


+  P2[(B-i-l)j  ij  {l-P**(ij)}  "(zjjdzj-zj  |l-P**(2j)}  ]. 


Upon  substituting  the  above  expression  of  K(zj^)  in  (6.29)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.28). 

Relation  6.11;  For  n>2, 

-  Clll} 

+  (»-»{ Vl,a:n  -  Vi-i}  -  •  («•«) 

Proof.  For  n>2,  let  us  consider 


**(1)  -  w(7*  7*0  X 

®n-l:n  '  n-l:n  ^n:n' 


=  n(n-l) 


(Zj)  K(zj)  dz^,  (6.32) 
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where 


K(zi)  ”  I  ^  ^*2^  ^2  ‘ 


(6.33) 


Upon  using  (6.6)  in  (6.33)  and  integrating  by  parts,  we  get 


Pi 

K(Zj)  =  (2P-1)[  I  f**(*2)  d22  -  *1  J 


Pi  2 

-  (P-Q)[2j  *2  f"(*2)  -  "1  {l-P  f^))  ] 


+  Pj  (J’l  -  *l)- 

Upon  substituting  the  above  expression  of  K(*i)  in  (6.32)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.31). 


1:  For  n>2. 


‘4!3:b+1  =  “si+l  5(P^[“2:i^'  '  “'^2(“li-l  " 

-  -  “1*2:4]  ■ 
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Proof.  For  n>2,  let  us  start  with 


*♦(1)  -W7*^  7*  \ 

“2:2  ”^^n:2^2:2^ 


—  2J  Z2  f  (zj)  K(z2)  dz2, 


where 


K(zJ 


(zj)  dzj 


(6.35) 


(6.36) 


Upon  using  (6.5)  in  (6.36)  and  integrating  by  parts,  we  get 


*2 

K(z2)  =  Q2(22“Qi)  (1-2Q){22  ^**(*2^  "  ^**(*l)  ‘^l} 


-  (P-Q)[z2  {^**(*2^}  -  2^  *1  F*%i)  f**(zj)  dzj  . 


Upon  substituting  the  above  expression  of  K(z2)  in  (6.35)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.34). 

Relation  6.13:  For  2<i<n-l, 


“i+14+2:n+l  "  “i+Ln+1  (I+l 


[1  i+l:n  i  (  i:n-l  i-l,i:n-lj 


(6.37) 
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Proof.  For  2<i<n-l,  let  us  consider 


*♦(1)  -F(7*0  7*  \ 


where 


(6.38) 


r2/-  4,4,  ii-1 

I  |f  (z^)|  f  (z^)  dz^  .  (6.39) 


Upon  using  (6.5)  in  (6.39)  and  integrating  by  parts,  we  get 


Upon  substituting  the  above  expression  of  K(z2)  in  (6.38)  and  simplifying  the  resulting 
equation,  we  derive  the  recurrence  relation  in  (6.37). 
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In  particular,  by  setting  i=n-l  in  (6.37),  we  obtain  the  recurrence  relation 


.  n+1  r  1  *♦(!)  **  1 

®n,n+l:n+l  ®n+l:n+l  n(P-Q)[n^  “n:n  n— l[“n-l:n-l  “n-2,n-l:n-lj 


(6.40) 


** 


**2 


It  should  be  mentioned  here  that  by  starting  with  the  result  that  2’2~^1'1 

(Govindarajulu,  1963;  Joshi,  1971),  one  may  employ  Relations  6.8  -  6.13  to  compute  all 

the  immediate  upper-diagonal  product  moments,  viz.,  i+i-n  ^  simple 

recursive  way  for  all  sample  sizes.  As  mentioned  earlier  in  Section  4,  this  is  sufficient  for 

*♦ 

the  evaluation  of  all  the  product  moments  as  the  remaining  product  moments,  viz.,  a. 
for  l<i<j<n  and  j-i>2,  may  be  determined  by  using  the  recurrence  relation  in  (4.13). 
However,  for  the  sake  of  completeness  we  present  here  some  more  recurrence  relations 
satisfied  by  the  general  product  moments.  These  results  may  be  established  by  following 
exactly  the  same  steps  as  used  in  proving  Relations  6.8  -  6.13. 


Relation  6.14:  For  l<i<n-2. 


** 


*i  ,n:n+l 


=  “i,n-l:n+l  +  '  “i.n-lm-l} 


Relation  6.15:  For  l<i<j<n-l  and  j-i>2, 


** 


** 


,  j:n+l  ®i ,  j-l:n+l 


n+1  “^2  1.**  \ 

(n-j+2){P-y)[H:Tfn“i , j:n-l  "  ®i  j 


(6.42) 
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Relation  6.16:  For  3<j<ii, 


“2,j+l:n+l  -  ^.j+lm+l  +  ‘ 


Relation  6.17:  For  2<i<j<n  and  j-4>2, 


(6.43) 


** 

^i+l,j+l:n+l 


*  * 

1 

i+2,j+l:n+l 

+  (i+l)(F-gjL^ 

“jrn  -r 

*♦  1 

**  11 

‘i-l,j-l:n-lj  ■ 

■(l~2Q)|“i+lJ:n 

”  “i ,  j:n|j  ’ 

(6.44) 


7.  DETAILS  OF  AVAILABLE  TABLES 

We  list  below  the  tables  that  are  currently  available  on  order  statistics  and  their 
moments. 

a.  Table  of  100a%  points  (for  o=0.50,  0.75,  0.90,  0.95,  0.975,  0.990)  for  all  order 
statistics  for  sample  sizes  up  to  10  and  for  extreme  and  central  order  statistics  for 
sample  sizes  from  11  to  25  has  been  given  by  Gupta  and  Shah  (1965); 

b.  Table  of  probability  integrals  of  the  sample  range  for  n=2  and  3  evaluated  at 
w=:0.20(0.20)1.00(0.50)4.00  has  been  given  by  Gupta  and  Shah  (1965); 
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c.  Table  of  means  and  standard  deviations  of  order  statistics  for  sample  sizes  up  to  10 
has  been  given  by  Birnbaum  and  Dudman  (1963).  Table  of  covariances  of  order 
statistics  for  sample  sizes  up  to  10  has  been  given  by  Shah  (1966).  These  two 
tables  have  been  extended  by  Gupta,  Qureishi  and  Shah  (1967)  for  sample  sizes  up 
to  25.  Recently,  Balakrishnan  and  Malik  (1990)  have  prepared  tables  of  means, 
variances  and  covariances  for  sample  sizes  50  and  less  in  which  the  values  are 
reported  to  ten  decimal  places; 

d.  By  using  the  results  presented  in  Section  6,  Balakrishnan  and  Joshi  (1983  b)  have 
given  tables  of  means,  variances  and  covariances  for  the  symmetrically  truncated 
logistic  distribution  (with  Q=l-P=0.01,0.05(0.05)0.20)  for  sample  sizes  up  to  10; 

e.  Balakrishnan  (1985)  has  handled  the  half  logistic  distribution  (case  when  and 
P=l)  and  has  presented  tables  of  means,  variances  and  covariances  for  sample  sizes 
up  to  15.  B:e  has  also  given  tables  of  100a%  points  (for  a=0.01,  0.05,  0.10(0.10) 
0.90,  0.95,  0.99)  for  extreme  order  statistics  for  sample  sizes  up  to  15.  In  addition, 
he  has  presented  a  table  of  modes  of  all  order  statistics  for  sample  sizes  up  to  15. 

8.  PLACKETT’S  APPROXIMATION 


David  and  Johnson  (1954)  and  Clark  and  Williams  (1958)  have  developed  some 
series  approximations  for  moments  of  order  statistics  from  an  arbitrary  continuous 
distribution.  These  have  been  developed  by  applying  the  probability  integral 
transformation  and  then  using  the  known  moments  of  order  statistics  from  the  uniform 
distribution.  Plackett  (1958),  instead,  has  used  the  logit  transformation  which  transforms 
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an  order  statistics  from  an  arbitrary  continuous  distribution  into  the  order  statistic 

-2 

from  the  logistic  L(0,^  )  distribution  to  develop  some  series  approximations  for  the 
moments  of  in  terms  of  the  moments  of  logistic  order  statistics  Zj.^. 

We  have  already  seen  in  Section  3  that  the  moments  and  the  cumulants  of  the 
logistic  order  statistics  Z|.^  are  all  available  in  explicit  form.  Now,  by  realizing  that  the 
logit  transformation 


Z  =  Ai 


iT=F^l 


(8.1) 


transforms  the  order  statistic  from  an  arbitrary  continuous  distribution  with  cdf 
Frp(t)  into  the  logistic  order  statistic  Z.  and,  therefore,  expanding  T.  in  a  Taylor 

X  l«ll  1»I1 

series  about  the  point  E(Zj.j^)  =  we  derive 


+ .. 


(8.2) 


where,  for  n-i+l>i. 


as  derived  in  Section  3,  and  t^^^  is  the  value  of  the  derivative  of  t  with  respect  to  Z  at 
Z  at  Z  =  l>y  taking  expectation  on  both  sides  of  (8.2)  and  upon  using 
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the  exact  and  explicit  expressions  for  the  cumulants  of  logistic  order  statistics  derived  in 
Section  3,  we  obtain  the  series  approximation 


e[t,J  =  .(0)  +  1  .P)  ^  ,  1  .(3)  .-(3)  ,  ^  .(4){,-(4),3[.-(3)]2}. 

(8.3) 


The  derivatives  appearing  as  coefficients  in  the  approximation  in  (8.3)  are  easy  to 
obtain  as  in  the  case  of  approximations  due  to  David  and  Johnson  (1954)  and  Clark  and 
Williams  (1958).  For  example,  for  the  standard  normal  distribution  with  probability 
density  function  <^(t)  and  cumulative  distribution  function  $(t),  we  have 


i:n 

t(l)  «  $  (1^)/(|», 

t(2)  ^  t(l)  I  ^(1)  _  (2$  _  1) 


t(3)  _  (  ^(1))  ^  2tt(^^  t^^)  +  t(^)  (1  -  2$)  -  2t(^4  (l-$), 


and 


t(^)  =  5(t(^))^  {2tt(^^  -  (2$  -  1)| 


+  2t(^)|tt(^^  -  2$  (1-#)|  +  2t(^)  (2$  -  1)  $  (1  -  $). 
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The  above  given  derivatives  are  all  bounded.  As  pointed  out  by  Blom  (1958),  suppose  we 
include  the  first  j-1  terms  in  the  series  expansion  for  E(T.,^)  obtained  from  (8.2),  then 
the  absolute  value  of  the  remainder  after  j-1  terms  is  at  most 


Since  ElZ.  „ 
'  i:n 


j  max  I  .0)  I 

'‘i'B  is  known  and  also  that 


(8.4) 


{  ®|Zi:n  -  ^  <  {  E|Z..„  -  p, 


we  will  be  able  to  present  bounds  to  E(T.  )  for  all  values  of  j. 
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